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ABSTRACT.  Guided  wave  propagation  in  a  homogeneous  plate  is  known  to  be  sensitive  to  both 
temperature  changes  and  applied  stress  variations.  Here  we  consider  the  inverse  problem  of 
recovering  biaxial  stresses  from  measured  changes  in  phase  velocity  at  multiple  propagation 
directions  using  a  single  mode  at  a  specific  frequency.  These  changes  depend  upon  both  the 
magnitude  and  orientation  of  the  principle  stresses.  Although  there  is  no  closed  form  solution,  prior 
results  indicate  that  phase  velocity  changes  exhibit  a  sinusoidal  angular  dependence.  Here  it  is 
shown  that  all  sinusoidal  coefficients  can  be  estimated  from  a  single  uniaxial  loading  experiment. 
The  general  biaxial  inverse  problem  can  thus  be  solved  by  fitting  an  appropriate  sinusoid  to 
measured  phase  velocity  versus  propagation  angle,  and  relating  the  coefficients  to  the  unknown 
stresses.  The  phase  velocity  data  are  obtained  from  direct  arrivals  between  guided  wave  transducers 
whose  direct  paths  of  propagation  are  oriented  at  different  angles.  This  method  is  applied  to 
experimental  sparse  array  data  recorded  during  a  fatigue  test,  and  the  additional  complication  of  the 
resulting  fatigue  cracks  interfering  with  some  of  the  direct  arrivals  is  addressed  via  proper  selection 
of  transducer  pairs.  Results  show  that  applied  stresses  can  be  successfully  recovered  from  the 
measured  changes  in  guided  wave  signals. 

Key  words:  Acoustoelasticity,  Lamb  waves,  Load  estimation 
PACS:  43.20.Bi,  43.35.Cg,  46.40.Cd 


INTRODUCTION 

Guided  waves  such  as  Lamb  waves  play  a  significant  role  in  nondestructive 
evaluation  (NDE)  and  structural  health  monitoring  (SHM)  techniques,  which  generally 
require  or  assume  insensitivity  to  varying  environmental  and  operational  conditions. 
However,  guided  ultrasonic  waves  are  well-known  to  have  unavoidable  sensitivity  to 
environmental  changes  such  as  temperature,  surface  wetting  and  applied  loads.  However, 
applied  loads  can  also  open  tightly  closed  fatigue  cracks,  so  it  is  of  interest  to  know  the 
current  loading  state.  This  paper  describes  an  inverse  method  to  estimate  the  stress  tensor 
as  a  homogeneous  biaxial  load  is  applied,  which  is  based  on  the  forward  problem  of 
calculating  dispersion  curves  for  acoustoelastic  Lamb  waves  [1].  The  principal  stress 
components  and  orientation  are  estimated  from  a  sinusoidal  fit  of  ultrasonic  data  collected 
from  the  same  spatially  distributed  array  that  is  being  used  to  detect  and  characterize 
damage.  The  proposed  method  is  experimentally  validated  using  fatigue  test  data  acquired 
on  an  aluminum  plate  from  an  array  of  spatially  distributed  piezoelectric.  To  minimize 
scattering  effects  caused  by  growing  fatigue  cracks,  a  subset  of  all  available  transducer 
pairs  is  selected  and  used  to  obtain  a  better  sinusoid  fit. 
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BACKGROUND 


Consider  a  homogenous,  isotropic  aluminum  plate  with  thickness  d  and  infinite  in 
extent  as  shown  in  Figure  1.  Biaxial  stresses  an  and  an  are  applied  along  the  x\  and  xi  axes 
of  a  rectilinear  coordinate  system  x,  =  (xi,  xo,  X3),  which  are  rotated  by  an  angle  a  from  a 
measurement  coordinate  system  indicated  by  ( xd,  X2',  X3').  Assume  that  ultrasonic  guided 
waves  are  propagating  along  a  direction  that  makes  an  angle  Q  with  respect  to  the  x\  axis 
and  0'  with  respect  to  the  xi'axis. 

As  per  previous  work,  the  theory  for  acoustoelastic  Lamb  wave  propagation  has 
been  developed  [1]  for  biaxial  loads  and  an  arbitrary  direction  of  propagation.  Results  show 
that  isotropic  dispersion  curves  for  a  stress-free  plate  become  anisotropic  with  phase 
velocities  becoming  angle  and  stress  dependent  for  a  specified  mode  and  frequency. 
Changes  of  phase  velocity  with  load  can  be  accurately  approximated  as  a  sinusoidal 
function  with  respect  to  the  propagation  angles  as  shown  in  Figure  2  for  the  fundamental 
symmetric  (So)  mode. 


FIGURE  1.  Geometry  for  Lamb  wave  propagation  under  applied  biaxial  stresses. 


(a) 


(b) 


FIGURE  2.  (a)  Dispersion  curves  for  the  S0  mode  when  cr,  ,  =  30  MPa,  a2 2  =  60  MPa  for  different  propagation 
angles,  (b)  Changes  of  phase  velocity  for  the  So  mode  when  an  =  30  MPa,  a22  =  60  MPa  at  400  kHz  as  a 
function  of  propagation  angle. 
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LOAD  ESTIMATION  PROCEDURE 


Based  on  the  theory  for  acoustoelastic  guided  wave  under  a  uniaxial  load,  the  phase 
velocity  change  A cp  for  a  given  frequency  and  Lamb  mode  has  the  following  form  [2] : 

Ac  I  =  <7n(X  cos2  0  +  K2  sin2  0)  (1) 

1  l<722=0 

A cp  |  o  =  (J22  (K3  cos2  0  +  K4  sin2  0)  (2) 

In  these  equations  an  and  <722  are  uniaxial  applied  stresses  in  the  x\  and  X2  directions, 
respectively,  6  is  the  direction  of  Lamb  wave  propagation  in  the  principal  (unprimed) 
coordinate  system,  and  K\,  K2,  K3  and  K4  are  the  four  acoustoelastic  constants  for  the 
particular  frequency,  mode  and  loading  direction.  These  equations  are  the  same  form  as  for 
non-dispersive  bulk  and  Rayleigh  waves,  although  the  acoustoelastic  constants  for  Lamb 
waves  are  frequency  and  mode-dependent  even  for  homogeneous  and  isotropic  media  [3]. 

A  combined  equation  for  an  arbitrary  biaxial  load  can  be  deduced  by  first  noting 
that  the  four  acoustoelastic  constants  can  be  reduced  to  two  (i.e.,  K\  =  K4,  K2  =  ^3)  because 
of  symmetry.  Next,  we  assume  a  combined  sinusoidal  function  to  describe  the  changes  of 
phase  velocity  as  a  linear  combination  of  the  two  uniaxial  loading  cases,  which  can  be 
expressed  as  follows: 


A cp(0)  =  (A'1(j11  +  K2<j22) cos2  6  +  (K2<ju  +  Kxg22) sin2  9 .  (3) 

As  written  above,  changes  in  phase  velocity  are  expressed  in  the  principal  axis  system, 
which  may  not  coincide  with  the  measurement  system,  which  is  rotated  by  a  from  the 
principal  axis  system.  After  some  algebra,  Eq.  (3)  can  be  rewritten  as  follows: 

A  cp(9')  =  {Klan  +  K2a22)cos\0'+  a)  +  (K2crn  +  Klcr22)sm2(8’+  a)  (4) 

=  a0  +a1cos(20')  +  a2  sin(20'). 

In  this  equation, 

ao  =\(KX  +  K2){an  +  (J22),  flj  =\(K{ -K2)(an -cr22) cos(2or),  and  (5) 

a2  =  -\(KX  -  K2)(on  -  G22)sm(2a) 

Numerical  simulations  were  performed  to  validate  Eq.  (4).  First,  phase  velocity 
changes  for  the  So  mode  at  400  kHz  were  simulated  for  multiple  uniaxial  stress  conditions 
for  an  =  0  and  cr22  varying  from  0  MPa  to  100  MPa  in  steps  of  10  MPa;  the  propagation 
angle  6  varied  from  0  degree  to  90  degrees  in  steps  of  5  degrees.  Second,  K\  and  K2  were 
estimated  by  least-squares  using  Eq.  (2)  for  the  multiple  known  uniaxial  cases;  their  values 
were  consistent  for  all  cases  considered.  Third,  theoretical  phase  velocity  changes  were 
calculated  as  a  function  of  propagation  angle  for  multiple  cases  of  a,  an  and  cr22  as 
described  in  [1],  Finally,  constants  K\  and  Ki  determined  from  the  first  step  were  used  in 
Eq.  (4)  to  estimate  the  changes  of  phase  velocity  for  the  multiple  cases  of  a,  an  and  cr22- 
The  results  of  Eq.  (4)  were  then  compared  to  calculated  curves  and  were  found  to  be  in 
excellent  agreement  for  all  cases.  Typical  results  are  shown  in  Figure  3  for  a  =  30  degrees, 
an  =  30  MPa  and  aoi  =  60  MPa.  Therefore,  the  assumed  sinusoidal  dependence  and  linear 
combination  of  uniaxial  loading  cases  for  biaxial  loading  is  numerically  confirmed. 
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(a)  (b) 

FIGURE  3.  Phase  velocity  changes  for  the  S0  mode  at  400  kHz  versus  propagation  angle,  (a)  crn  =  30  MPa, 
cr22  =  60  MPa,  and  a  =  30  degrees,  (b)  a\  ,  =  30  MPa,  a 22  =  60  MPa,  and  a  =  60  degrees. 


The  approach  to  the  inverse  problem  can  be  readily  seen  by  considering  Eq.  (4),  which 
described  the  expected  sinusoidal  form  of  the  A cp  vs.  6 '  data.  Once  K\  and  Ki  are  obtained, 
the  constants  ao,  ti\  and  a  2  can  be  determined  via  least-squares.  Finally,  the  unknown 
biaxial  loads  on  and  an  along  with  the  orientation  angle  a  can  be  expressed  in  terms  of  ao, 
a  1  and  ao  as  follows: 


°ii  = 


a0cos(2a){Kl  -  K2)  +  al{Kl  +  K2) 


cos(2  CC)(K[ 


■K) 


and  a  =  —  arctan 
2 


1  J 


aQcos(2a)(K2  -  Kl)  +  al(Kl  +K2) 
cos(2  a)(K2-Kf) 


(5) 


This  procedure  is  used  in  the  subsequent  sections  to  estimate  stresses  from  measurements 
of  phase  velocity  changes  with  load. 

EXPERIMENTAL  VALIDATION 

A  fatigue  test  was  performed  with  an  array  of  six  surface-bonded  PZT  transducers 
on  a  6061  aluminum  plate  as  shown  in  Figure  4.  The  specimen  was  fatigued  using  a 
sinusoidal  tension-tension  profile  from  16.5  MPa  to  165  MPa  at  3  Hz.  Ultrasonic  signals 
from  the  15  unique  transducer  pairs  were  recorded  for  uniaxial  loads  ranging  from  0  to 
115  MPa  in  steps  of  1 1.5  MPa  for  each  data  set.  A  total  of  fourteen  data  sets  were  recorded, 
where  each  data  set  contains  11  static  loading  measurements.  Additional  information, 
including  the  growth  of  fatigue  cracks,  is  summarized  in  [4]. 

Guided  waves  were  generated  by  a  broadband  chirp  excitation.  Then,  the  measured 
signals  were  filtered  using  a  7  cycle,  Hanning  windowed,  400  kHz  tone  burst  signal  [5]. 
In  this  paper,  the  So  Famb  wave  mode  was  chosen  for  analysis  because  it  has  clear  first 
arrivals  from  all  the  transducer  pairs,  which  is  more  convenient  than  the  slower  Ao  mode 
for  accurately  extracting  phase  velocity  changes. 
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FIGURE  4.  Drawing  of  the  specimen  and  transducer  geometry  (not  to  scale). 


FIGURE  5.  (a)  First  arrivals  of  transducer  pair  #2-5,  data  set  #1,  for  the  11  uniaxial  loading  conditions 
(baseline,  no  crack),  (b)  Zero  crossing  times  with  respect  to  loads  for  transducer  pair  #2-5,  data  set  #1 . 

Changes  in  received  signals  from  each  transducer  pair  caused  by  uniaxial  loading  were  first 
investigated  by  examining  zero  crossing  times  within  the  first  arrival  of  each  signal.  For 
example,  Figure  5(a)  shows  received  signals  from  transducer  pair  #2-5  and  data  set  #1  (no 
cracks)  for  the  11  uniaxial  loads.  Small  time  shifts  of  the  first  arrival  with  each  load 
increment  can  be  seen,  and  a  linear  relationship  between  zero  crossing  time  and  applied 
load  is  observed  as  shown  in  Figure  5(b).  Changes  in  phase  velocity  for  the  11  loading 
conditions  and  15  transducer  pairs,  which  refer  to  different  propagation  angles,  were 
extracted  from  first  arrivals  and  zero  crossing  times.  Constants  K\  and  Ki  were  estimated 
from  all  transducer  pairs  and  known  loading  conditions  of  data  set  #1  via  Eq.  (2). 

Loading  conditions  were  assumed  to  be  unknown  for  all  other  data  sets.  Constants 
ao,  ci\,  and  go  were  estimated  via  a  sinusoidal  least  squares  fit  of  phase  velocity  changes  vs. 
propagation  angle  using  Eq.  (4),  and  an,  022  and  a  were  calculated  from  a0,  ci\,  and  c/2 using 
Eq.  (5).  As  an  example,  Figure  6  shows  the  sinusoidal  fit  for  changes  in  phase  velocity  with 
respect  to  the  propagation  angle  using  all  the  transducer  pairs  of  data  set  #2  for  two 
different  loading  conditions.  Table  1  shows  that  the  recovered  loads  and  orientation  angle 
are  in  good  agreement  with  actual  values. 
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(a)  (b) 

FIGURE  6.  Example  experimental  data  and  sinusoidal  fit  for  data  set  #2.  (a)  an  =  0  MPa,  a2 2  =  57.5  MPa. 
(b)  (Tu  =  0  MPa,  cr22=  115  MPa. 


TABLE  1.  Estimated  stresses  and  angle  for  two  loads  from  data  set  #2 


(Tu  (MPa) 

022  (MPa) 

a  (degree) 

0ii  (MPa) 

022  (MPa) 

a  (degree) 

Actual 

0 

57.5 

0 

0 

115 

0 

Estimated 

-1.24 

55.65 

-0.15 

-2.20 

112.07 

-0.03 

However,  as  cracks  appear  and  grow  in  the  later  data  sets,  the  ultrasonic  signals 
obtained  from  the  transducer  pairs  whose  direct  paths  travel  through  cracks  become  much 
more  complicated.  As  an  example,  Figure  7(a)  shows  both  nonlinear  time  shifts  of  the  first 
arrivals  and  significant  amplitude  decreases  from  data  set  #10  (two  cracks)  because  the 
fatigue  cracks  interfered  with  the  direct  ultrasonic  arrivals.  Figure  7(b)  shows  that  when 
the  two  cracks  open,  the  zero  crossing  times  are  no  longer  linear  with  respect  to  the  applied 
load,  which  will  clearly  affect  the  accuracy  of  the  computed  changes  in  phase  velocity  with 
load. 


FIGURE  7.  (a)  First  arrivals  of  11  uniaxial  loading  conditions  for  transducer  pair  #2-5,  data  set  #10  (two 
cracks),  (b)  Zero  crossing  times  with  respect  to  loads  for  pair  #2-5,  data  set  #14. 
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The  complication  of  cracks  interfering  with  the  direct  arrivals  is  mitigated  by 
excluding  some  of  the  signals  from  the  sinusoidal  fit.  Referring  to  Figure  4,  it  can  be  seen 
that  the  direct  arrivals  of  three  transducer  pairs  travel  directly  through  the  cracked  region: 
#1-4,  #2-5  and  #3-6.  If  these  pairs  are  excluded,  the  effects  of  cracks  are  minimized. 
Figure  8  shows  the  sinusoidal  fit  for  changes  in  phase  velocity  with  respect  to  the 
propagation  angle  from  two  loading  conditions  using  both  the  12  selected  pairs  and  all  15 
pairs. 

Finally,  an  and  a  in  along  with  the  angle  a  were  calculated  from  a(h  a\,  and  cii  from 
Eq.  (5).  Figure  9(a)  shows  the  estimated  uniaxial  loads  and  angle  when  actual  <j\  \  =  0  MPa 
and  <722=  57.5  MPa,  and  Figure  9(b)  shows  the  results  for  on  =  0  MPa  and  022  =  115  MPa. 
Results  for  both  the  selected  pairs  and  all  the  pairs  are  shown  for  comparison,  where  it  can 
be  seen  that  excluding  some  of  the  transducer  pairs  considerably  improved  results  for  the 
later  data  sets. 


(a) 


(b) 


FIGURE  8.  Experimental  data  and  sinusoidal  fit  of  phase  velocity  changes  for  data  set  #10.  (a)  on  =  0  MPa, 
o-22  =  57.5  MPa,  (b)  an=0  MPa,  a22  =  1 15  MPa. 


Data  Set  Number 


(a) 


Data  Set  Number 


(b) 


FIGURE  9.  Estimated  stresses  and  orientation  angles  for  all  data  sets,  (a)  an  =  0  MPa,  er22  =  57.5  MPa, 
a  =  0  degrees,  (b)  an  =  0  MPa,  er22  =115  MPa,  a  =  0  degrees. 


7 

Approved  for  public  release;  distribution  unlimited. 


CONCLUSIONS 


This  paper  shows  a  load  estimation  strategy  based  on  the  assumption  that  the  Lamb 
wave  acoustoelastic  response  from  a  biaxial  loading  case  can  be  decomposed  into  that  from 
the  two  orthogonal  uniaxial  loading  cases.  This  assumption  is  verified  numerically,  and  it 
is  further  shown  that  only  two  acoustoelastic  constants  are  needed  to  describe  the  general 
angle  and  stress  dependence  of  a  specific  Lamb  wave  mode  and  frequency  to  a 
homogeneous  biaxial  load.  Unknown  applied  stresses  and  direction  can  be  estimated  using 
the  two  acoustoelastic  constants  by  finding  the  best  sinusoidal  fit  of  phase  velocity  changes 
versus  propagation  direction.  The  efficacy  of  the  proposed  strategy  is  verified  by  a  fatigue 
test  having  multiple  ultrasonic  measurements  from  various  uniaxial  static  loading 
conditions.  A  reduced  set  of  transducer  pairs  was  selected  to  improve  the  accuracy  of  the 
estimation  results  by  minimizing  the  interference  of  the  direct  ultrasonic  waves  with  fatigue 
cracks  as  they  open  under  applied  loads. 
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